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1. INTRODUCTION 

Variability of the spectral energy distribution (SED) of active galactic 
nuclei (AGNs) can provide clues for understanding both the main emis- 
sion processes and the origin of their variations. So far, multi-wavelength 
monitoring has been possible for a limited number of bright objects, 
thanks to large international cooperations (see Ulrich, Maraschi, & Urry 
1997, and refs. therein). The "ensemble" analysis of the light curves 
was obtained in the past only from single-band observations of statisti- 
cal samples of QSOs (Giallongo, Trevese, & Vagnetti 1991; Hook et al. 
1994; Cristiani et al. 1996; Di Clemente et al. 1996). Cutri et al. (1985), 
Kinney et al. (1991), and Edelson et al. (1990), found a hardening of 
the spectrum in the bright phase in small QSO samples. A positive 
correlation of the spectral index a {f v oc u a ) with brightness variations 
A log f u has been found by Trevese et al. (2000). The recent publication, 
by Giveon et al. 1999, of two-band monitoring of 42 PG QSOs, allows 
a statistical study of the spectral variability of the entire set of light 
curves (Trevese & Vagnetti 2001). They find a positive correlation of 
the B — R color with the brightness variations 5B. These results confirm 
the interpretation of the variability-redshift (v-z) correlation suggested 
by Giallongo, Trevese, & Vagnetti (1991) and Di Clemente et al. (1996). 
An independent analysis of the same light curves, performed by Cid 
Fernades et al. (2000) implies the existence of an underlying spectral 
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component, redder than the variable one, which could be identified ei- 
ther with the host galaxy, considered by Romano & Peterson (1998), or 
with a non-flaring part of the QSO spectrum. We evaluate the effect of 
the host galaxy constant SED on the variable QSO+host SED, through 
numerical simulation based on templates of the QSO and host SEDs 
provided by Elvis et al. 1994. We also compare the results with the 
spectral variations deduced from Giveon et al. (1999), concluding that 
the nuclear spectral variability is dominant. 

2. DISCUSSION AND CONCLUSIONS 

From an atlas of 47 normal QSO continuum spectra, Elvis et al. 
(1994) derive the average QSO template spectrum corrected for the con- 
tribution of a template host galaxy spectrum. On the basis of these 
templates we produce a synthetic (QSO+host) SED. We characterize 
the relative contribution of the nucleus and the host by the parameter: 
T] = log(L^ / L g H ) , where L^j- and L 9 H are the total H band luminosities 
of the QSO and the host galaxy respectively. To test (disprove) the 
hypothesis that the spectral energy distribution of the active nucleus 
maintains its shape while changing brightness, we represent variability 
by small changes of Arj(t) around each adopted r] value. For each syn- 
thetic spectrum we compute the instantaneous spectral slope as a func- 
tion of frequency and time: a(u,t) = SlogL^/SlogzA We characterize 
the continuum SED variability by the (3(v) = da{y) / d\ogL u parameter 
representing the variation of the local spectral slope per unit change of 
log L u . 

Variability of QSO spectra has been analyzed on the basis of the data 
of Giveon et al. (1999), consisting of the light curves of a sample of 
42 nearby and bright QSOs (z < 0A,B < 16 mag), belonging to the 
Palomar-Green (PG) sample, monitored for 7 years with median sam- 
pling interval of 39 days, in B and R, with the 1 m telescope of the 
Wise Observatory. For each time interval r between two observations, it 
is possible to define (3(t) = [a(t + r) — a(t)] / [log Lqj(t + r) — log Lp-(t)], 
v = \fvBVR- We consider, for each object, the spectral variability param- 
eter [3 m defined as the average of (3(t) values for all r < 1000 days and 
the average value of the spectral slope a m . Both these quantities can be 
compared with the values (/3,a) obtained from the synthetic spectrum. 
In Figure 1 the continuous curve represents (5 versus a for r/ ranging 
from -3 (dominant host galaxy, in B and R) to 3 (~ pure QSO). To 
evaluate the effect of redshift we also show the dashed curve computed 
for z = 0.4, the maximum redshift of the Giveon et al. (1999) sample. 
On average, the distribution of points is clearly shifted towards top right 
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of Figure 1. This means that the relative weight of nuclear and stellar 
component cannot account for the observed spectral variability, which 
therefore must be intrinsic of the active nucleus. 
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Figure 1 The observed /3 
versus a values for PG 
QSOs. The two curves, 
(z = continuous, z=0.4 
dashed), represent /3(a) of 
the composite spectrum, 
for -3 < 77 < +3 . 



We have computed synthetic QSO+host SEDs on the basis of tem- 
plates of the nuclear and the star light spectra. We have computed the 
spectral slope a and the spectral variability (3 for the sample of PG 
QSOs monitored in the B and R bands by the Wise Observatory group. 
A comparison of the observed distribution of points in the (3 — a plane 
with the synthetic j3{a) allows to disprove the hypothesis that the ob- 
served variations of the QSO spectral shape are solely due to a change of 
the relative weights of the nuclear and stellar component. In the frame- 
work of the analysis of Cid Fernandes et al. (2000), our result implies 
that the constant spectral component should be identified with a non 
flaring part of the accretion disk. 
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